Mon. Not. R. Astron. Soc. 000, [T]-?? (2008) 



Printed 14 November 2011 



(MN XAT^ style file v2.2) 



Low frequency radio spectrum of LS 5039 during 
periastron and apastron passages 

Subir Bhattacharyya^* Sagar Godambe^ Nilay Bhatt^ Abhas Mitra^ Manojendu Choudhury' 

^ Astrophysical Sciences Division, Bhabha Atomic Research Centre, Mumbai 400085, India 
^ Homi Bhabha Centre for Science Education, TIFR, Mumbai 400094, India 

14 November 2011 

ABSTRACT 

We have recently studied LS5039, a gamma-ray binary, with Giant Meterwave 
Radio Telescope (GMRT) during its periastron and apastron passage. The results 
presented here show that the spectra are inverted at the low frequency and the flux 
densities do not differ significantly for two different orbital phases. Assuming that 
the free-free absorption of radio in stellar wind environment is responsible for the 
optically thick radio emission we calculated the free-free absorption optical depth and 
constrained the height of the radio emitting region from the orbital plane. The height 
is found to be around 1.6 AU for a spherical stellar wind geometry. This estimate may 
change if the stellar wind is focussed or the radio absorption is due to synchrotron 
self-absorption. 
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1 INTRODUCTION 

The major unresolved issues regarding the gamma-ray bi- 
nary LS 5039 are (i) the nature of the compact object is not 
known, (ii) till date there is no convincing evidence in favour 
of the existence of any accretion disc in the source, and (Hi) 
the nature of th e stellar wind is not clearly understood. 
iLiu et all l|2006l ') showed LS 5039 to be a high mass X-ray bi- 
nary with a ON6.5V( f) ty pe primary star of ma ss 22.9 Mq. 
iMcSwain et~all l|2004 l and lCasares et al.l l|2005f ) studied the 
source i n UV and optic a l. Th rough UV and optical spec- 
troscopy [McS wam et all l|2004 ) found P Cygni line profiles 
in the UV Nv A 1204 an d CIV A 1550 l ines in dicating the 
presence of strong wind. ICasares et al.l (|2005l ) studied the 
optical H Balmer and Hel and Hell lines and determined 

the orbital parameters. 

In a recent study carried out bv lSartv et al.l l|201ll ) using si- 
multaneous observations with the space based Canadian Mi- 
crovariability and Oscillation of Stars (MOST) and ground 
based 2.3 m optical telescope Australian National University 
(ANU), the stellar wind properties and the orbital parame- 
ters were studied extensivel y. Orbital parameter s were better 
constrained as compared to lCasares et al.l (|2005l ) . The eccen- 
tri city of the orb i t is fo und to be 0. 24 as comapred to 0.3 5 
bv lCasares et all (|2005l ) and 0.48 bv lMcSwain et all (|2004h . 
Since the compact object mass was found to be > 1.8Mq, 
we find it is difficult to conclude about the nature of the 
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compact object because with favourable equation of states 
the masses of neutron stars can easily approach such a value. 

Fro m their obser vations of H/? and Hel spectral lines 
iSarty et al.l l|201 j) revealed that equivalent width (EW) 
varies with the orbital phase and, particularly, EW is lower 
at inferior conjuction. This observation suggests that the 
stellar wind is possibly focussed towards the compact ob- 
ject instead of being spherically sym metric. A focuss e d win d 
is also observed in case of Cyg X-1 (|Gies fc BohonI l|l986t ): 
iMiUer et al.l l^2005^'). a high mass X-ray binary and it was 
modelled by ICies fc BoltonI (Il986l ). Even though the pres- 
ence of focussed stellar wind in LS 5039 needs further study 
and is to be modelled properly, but it possibly has direct 
implications on the results we present in this paper. 

In our previous work (|Godambe et al.l (|2008l ). Paper I), it 
was shown that the spectrum in the low frequency range is 
optically thick and there is a possible spectral turn over 
at around 1 GHz. The spectral index was found to be 
—0.749 ± 0.111. The slope of the optically thick spectrum 
differs from the ideal condition of synchrotron self-absorbed 
spectrum which shows a spectral slope of |. The opti- 
cally thick spectrum in radio can also be generated by free- 
free absorption of radiation in presence of the stellar wind 
from the hot and massive companion star. Considering the 
above two pos s ible sc enarios of absorption of radio spectrum 
iBosch- RamonI (|2009l ) constrained the physical conditions of 
the radio emitting region in LS5039 jet. In case of free- free 
absorption of radio wave in the stellar wind of the compan- 
ion the magnetic field in the emission region was estimated 
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Table 1. Log of GMRT observations 
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Frequency 


Phase 




(MHz) 


range 


54573 
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0.63 - 0.70 




605.2 
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Figure 1. Radio spectrum of LS 5039 as observed with GMRT 



to be in the range (3 — 30) x 10~^ G while the size and 
the location of the emitter will be (3 — 4.5) x 10^^ cm and 
4.5 X 10^^ cm respectively. If SSA was responsible for the 
optically thick radio spectrum then the magnetic was found 
to be in the range 10"^ — 1 G. The size of the emitting region 
would be ~ (4.5 — 10) x 10^'^ cm and location of the source 
would be at > (4.5 — 10) x 10^^ cm. Hence it is important to 
understand the actual process responsible for the absorption 
of radio emission. 

Here we present the results of radio observations of LS 5039 
taken with Giant Meterwave Radio Telescope in the 200 
- 1280 MHz frequency range. We study the low frequency 
spectrum during the periastron and apastron passage of the 
compact star and attempt to understand its possible conse- 
quences. The observation and the data analysis are discussed 
in the next section. The results are discussed in Section 3 
and finally we conclude the paper in Section 4. 



2 OBSERVATIONS AND DATA REDUCTION 

The GMRT consist of 30 steerable antennas of 45 m diam- 
eter in an approximate 'Y' shape similar to the VLA but 
with each antenna in a fixed position. Fourteen antennas are 
randomly placed within a central 1 km x 1 km square (the 
'Central Square') and the remaining antennas form the irreg- 
ular Y shape (six on each arm) over a total extent of about 
25 km. We refer the reader for details about the GMRT ar - 
ray to|http://gmrt.ncra.tifr.res.in and lSwarup et al.1 l|l99lf ). 



We observed LS 5039 at 1280 MHz on February 23, 2006 and 
simultaneously observed at 234 and 614 MHz on March 13, 
2006. 

The observations were made in standard fashion, with 
each source observation (30 min) interspersed with obser- 
vations of the phase calibrator (4 min) . The primary flux 
density calibrator was either 3 C 48 and 3C 286, w ith all flux 
densities being on the scale o fS aars et al.l j 19771 ). Either of 
these flux calibrator was observed for 20 min at the begin- 
ning and end of each observing session. 

The data recorded with the GMRT was converted to 
FITS format and analyzed with the Astronomical Image 
Processing System (AIPS) using standard procedures. Dur- 
ing the data reduction, we had to discard a fraction of the 
data affected by radio frequency interference and system 
malfunctions. After editing the data were calibrated and col- 
lapsed into a fewer channels. A self-calibration on the data 
was used to correct for the phase related errors and improve 
the image quality. The observation dates, the observed fre- 
quencies and the corresponding spectral phases are given in 
Tabled 



3 RESULTS AND DISCUSSION 

The spectra obtained for the two orbital phase ranges corre- 
sponding to the periastron and apastron passage are shown 
in Fig. [T] along wi th the data points from our previous obser- 
vations (Pa per I, Godambe et all (|2008|)) and archiva l data 
points from lMartf et al.1 (ll998l)andlRib6 et al.1 (|l999l ). The 
flux points obtained from 'Ribo et al.' ('1999') do not match 
with the points obtained by Marti ct al. (1 998). This is be- 
cause the observation by iRibo etld] (| 19991 ) with VLBI re- 
solved the source and we have used the flux values reported 
for the core. The overall flux from the core and the wings 
is of the same order as reported bv lMartf et al. I (|l998l ). It is 
evident that the spectra remain optically thick and the flux 
values at frequencies 605 and 1280 MHz are almost equal 
during the periastron (phase range : 0.65 - 0.72) and apas- 
tron (phase range : 0.99 - 0.06) passages. Because of poor 
quality of data at 234 MHz, we could not produce source 
image for both the observing spells. Compared to the flux 
values of our previous observations, the present flux values 
fall within the error bars of the previous observations. It is 
to be noted that the exposure time for our previous observa- 
tions was less compared to that of the present observations. 
The optically thick spectrum in radio frequencies can be pro- 
duced due to the absorption of radiation either by free-free 
absorption or by synchrotron-self absorption in the emitting 
region of the source. The free-free absorption of radio waves 
may occur in LS 5039 due to the presence of stellar wind fron 
the primary star. The relativistic jet where the radio emis- 
sion takes place by synchrotron process, moves through an 
environment of such a stellar wind. For the time being, let us 
assume that the wind from the companion star is spherically 
symmetric. If so, the f r ee-fre e absorption opacity is given by 
(Rvbicki fc LightmanI (| 19861 )) 



al^ = 0.018T' ^ Z^n^v' 



2- 

9ff 



(1) 



where T is the temperature of the stellar wind, n is the 
number density of the stellar wind and v is the frequency of 
radiation. 
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Figure 2. Schematic diagram of binary orbit and the jet 



For a binary orbit inclined at an angle i with respect to the 
sky plane, the jet will make an angle i with the ditection of 
observation. It is assumed that the jet is perpendicular to 
the orbital plane. If the semi- major axis of the orbit is a and 
the eccentricity is e then the equation of the orbit is given 
by 



(1 



(2) 



where 9 is the phase angle and 6 = 2n(f), ^ <^ ^ 1 is orbital 
phase. From Figure (2] using simple geometry, one can write 



O'C^ = + a'^e^ — 2pae cos t 



(3) 



If the radio emitting region is located at a height z from the 
base of the jet then 

= O'C^ + z'^ = + a'e^ — 2pae cos 9 + z'^ (4) 

In Figure El ZO'fC = (f - x), so 

— + ~ 2sx cos ('^ + (x ^ (5) 
Using equations ((2| and (U} equation ([Sjl can be written as 



1 — cos^ 



+ ae' - 2 



1 — cos^ 

2 , 2 



ae cos 9 + 



z -^^ X + 2sa:: sin(x — j) (6) 

Expressing all distances in the units of a, equation Q can 
be rewritten as 

f ^ = (9) + z'^ + x'^ + 2sx sin(x - i) (7) 

where f=-,z=-,x = - and 

K\9)= ^7^'2. +e'-2e/ ^~^\. ^os9 (8) 
1 — cos^ 9 V 1 — cos^ 9 

Using the geometry of Figure [2] one can finally write 

= + a^f + (9) 

where iF = z cos i — K{9) sini and SS = z sin i + K(9) cos i. 
For an aligned jet, where i = 0, 5" = 2 and ^ = K(ff). 
For a spherical wind the wind den sity at a distance r from 
star is given by (jPuls et al.l (|l996l )'). 



M 



M 



(10) 



Here M is the mass loss rate of the wind, is the radius of 



the star, Voa is the terminal velocity of the wind at far away 
distance. The free-free absorption optical depth is given by 



El 



dx. 



(11) 



As r > i?*, so with the approximation that (1 — 4f ) 
(1 + 2^), the optical depth is given by 



— sm I 2 tan 



2,:^3 



— — tan ^ i ^ 
2 



+ 



2R^ 



1 . 3 

+ - sm« 



tan-i ( :C 



— — sm I tan 



(12) 



For LS 5039, the parameter values are T = 37500K, 
M = 10"^ M© yr-\ a = 20 Rq, e = 0.24, = 9 5Rp ,, 
= 24° and Woo = 2440 km s"^ llMcSwain et al.1 (|2004l ): 



ICasares et"al] ()2005l l: [Mc"Swain et al.1 ()201lf )l. The variation 



of optical depth with orbital phase for three different lo- 
cations of the radio emitting regions are shown in Figure 
[3] for frequencies 1280 and 605 MHz. For the inner jet it is 
evident that the optical depth varies periodically. It is maxi- 
mum during the periastron passage (0 = 0) and is minimum 
during apastron passage (0 = 0.5). This is true for both 0° 
and 24° inclination angle of the jet vis-a-vis the inclination 
of orbital plane with respect to the sky plane. As the height 
of the radio emission zone increases the periodic variation 
of the optical depth diminishes and essentially the optical 
depth becomes almost independent of the orbital phase for 
2 « 17 (z = 17a = 1.6 AU). For z x 1.6 AU, the opti- 
cal depth for 1280 MHz, r//(!/ = 1280MHz) ^ 1 for the 
inclination angle 24°. As the spectral turnover is found to 
be around 1000 MHz, so r// fa 1 is the necessary condi- 
tion to constrain the emission height for the given physical 
condition. Since the optical depth at this height becomes 
independent of the orbital phase so the emitted radio flux 
densities will also be phase-independent. The approximately 
constant observed flux during periastron and apastron pas- 
sages reported here does not contradict the present absorp- 
tion model. For an aligned jet (i = 0), the estimated height 
of the radio emission zone will be little less than 1.6 AU. 
These estimates indeed assume that the basic emission pro- 
cess remains unchanged for all orbital phases. It is only the 
wind density changes for different orbital phases offering dif- 
ferent optical depth to the emitted radiation. 
However these estimations may change if (i) the stellar 
wind is focussed or (ii) the absortion is due to synchrotron 
self-absorption in the jet. In fact recent observation of He 
I and H/3 lines from LS 5039 by ISartv et al.l l|201ll ) indi- 
cates the p r esence of focussed wind in the binary system. 
Sart v et al.l (|201lf ) observed significant changes in the H/3 
and He I spectral lines with the orbital phase. They at- 
tributed the low equivalent width of the spectral line to the 
focussing of the stellar wind towards the compact object 
due to its strong gravity. iBosch-Ramon et al.l l|2007l ) anal- 



^ Similar estim ation of free-free optical depth is also done by 
IZdziarskil 1I2OIII') . 



Suhir Bhattacharyya Sugar Godambe Nilay Bhatt Abhas Mitra Manojendu Choudhury 
1280 MHz 605 MHz 



1.10 
1.06 
1.00 
0.95 
0.90 
0.85 
0.80 



— 1 — I — I — I — I — I — I — I — r 

z/a=17.0 



H — I — I — I — I — I — I — I — h- 

z/a = 5.0 



I I I I I I I I I 




-1 — I — I — I — I — I — r- 



I I I I I I I I I 



I I I I I I I I I 



170 
165 
160 
155 
150 
145 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



orbital phase ( 



Figure 3. Variation of free-free absorption optical depth with or- 
bital phase for 605 MHz and 1280 MHz for three different heights 
of the radio emitting region in the jet. Red curve: i = 24° and 
Green curve : j = 0°. Here the semi-major axis a is taken as 2QRq 



ysed the XMM-Newton data of LS 5039 during the perias- 
tron and apastron passage to study the effect of steUar wind 
on the X-ray absorption properties in this source. The spec- 
tral fitting of the data gave a value of equivalent hydrogen 
column density {Nh) which is co nsistent with the inte r stellar 
value, for all observation spells. iBosch-Ramon et al.l (|2007l ) 
used the photoelectric absorption model phabs. We reanal- 
ysed the same XMM-Newton data and fitted the spectrum 
with powerlaw and the improved X-ray absorption model 
tbabs using the spectral analysis tool, XSPEC. The X-ray 
absorption cross-section used in the model tbabs includes 
the cross-section for X-ray absorption by gas-phase ISM, by 
grain-phase ISM and by the molecules in the ISM. The fit- 
ted parameters for different observations are tabulated in 
Table (2] The value of Nh obtained in our fitting is compat- 
ible to the ISM value. This again indicates that the stellar 
wind does not contribute to the photoelectric absorption of 
the X-rays. This could be possible if the wind geometry is 
not spherical. Thus if the wind is focussed then the radio 
emission may take place at a lower height compared to the 
spherical wind scenario. So it is important to confirm the 
nature of stellar wind geometry for the companion star in 
LS 5039. In this context it is important to mention that Cyg 
X-1, a high mass X-ray binary, also exhibits focusse d ste l- 
lar wind fro m massive companion (Gics & Bolton ||1986D : 
iMiller et al.l |2005)). 

Another possibility of absorption of the radiation is due to 
synchrotron self- absorptio n of radio emi s sion i n jet. This 
scenario was discussed by i Bosch-RamonI (|2009l ) also. The 
absence of any orbital modulation of low frequency radio 
spectrum rule out this picture for th e inverted radi o spec - 
trum of LS 5039. The estimation by iBosch-RanionI (|2009l ') 
assumed an ideal | slope in the self-absorbed region of the 
spectrum. But observation shows that the spectral slope of 



optically thick portion of the spectrum deviates substatially 
from I (Paper I). It requires the modeling of the spectrum 
possibly generated in an inhomogeneous jet which is beyond 
the scope of this paper. To estimate the physical parame- 
ters of the jet it is important to model the multiwavelength 
data from the source. Such a model, incorporating the data 
from the low-frequency radio to gamma-rays, will be used 
to self-consistently produce the inverted radio spectrum, to 
reproduce the spectral turn over at around 1 GHz and to fit 
the broadband spectra, will be discussed in our next work. 



4 CONCLUSION 

We observed the gamma-ray binary LS 5039 in the low fre- 
quency radio by using Giant Meterwave Radio Telescope 
(GMRT) during the periastron and apastron passage of the 
source. Observations were carried out at frequencies 234, 605 
and 1280 MHz. The data at 234 MHz were discarded due 
to its poor quality. The analysed results show that the fiux 
values at 605 and 1280 MHz are consistent with our previous 
observations and they do not vary with the orbital phase. 
The spectrum is found to be inverted at low frequencies indi- 
cating an optically thick radio spectrum. Present data, when 
combined with the archival data, indicate a spectral turn 
over around IGHz (Figure 1). To understand the constancy 
of flux during apastron and periastron passages and to con- 
strain the location of radio emitting region we estimated the 
free-free absorption optical depth at radio frequencies due 
to the spherical stellar wind from the companion star in the 
binary system. From the free-free absorption optical depth 
estimates it is argued that the height of the radio emitting 
region along the jet is around 1.6 AU. 

It is to noted that the observed results presented here do 
not rule out the other possiblity of radio absorption by syn- 
chrotron self-absorption (SSA) process. If SSA is the respon- 
sible process then the above estimates may change. 
Also the estimation of the height of the radio emitting re- 
gion as described above assumes a spherical geometry of the 
stellar wind. Recent obser vations indicate th e presence of 
focussed wind in LS 5039 jSartv et all HoTJ)). If a model 
for focussed wind in this source is evolved and used to es- 
timate the free-free absorption optical depth then also the 
above estimates may change. Therefore to develop better 
understanding of radio emission and absorption processes in 
LS 5039 it is important to have more simultaneous multi- 
wavelength observations of the source. 
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Table 2. Spectral parameters for XMM-Newton observations for two different absorption models 



MJD 


ObsID 


Phase 


tbabs *powerlaw 


wabs *powerlaw 


photon index 


Nh X 10^2 cm-2 


photon index 


Nh X 10^2 cm-2 


54573 


0151160201 


0.53 


1.50 ± 0.03 


0.62 ±0.02 


1.49 ± 0.03 


0.63 ±0.03 


54670 


0151160301 


0.54 


1.48 ± 0.04 


0.66 ±0.03 


1.47 ±0.04 


0.67 ±0.03 


54684 


0202950201 


0.51 


1.53 ± 0.03 


0.69 ±0.02 


1.52 ±0.03 


0.71 ±0.02 


54688 


0202950301 


0.03 


1.63 ±0.04 


0.69 ±0.03 


1.62 ± 0.04 


0.71 ± 0.03 
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